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Abstract 

We apply the principles discussed in earlier papers to the construc- 
tion of discrete time quantum field theories. We discuss some of the 
issues to do with loss of Lorentz covariance and its recovery in the 
appropriate limit. We use the Schwinger action principle to find the 
discrete time free field commutators for scalar fields, which allows us to 
set up the reduction formalism for discrete time scattering processes. 
Then we derive the discrete time analogue of the Feynman rules for a 
scalar field with a cubic self interaction and give examples of discrete 
time scattering amplitude calculations. We find overall conservation of 
total linear momentum and overall conservation of total 6 parameters, 
which is the discrete time analogue of energy conservation and corre- 
sponds to the existence of a Logan invariant for the system. We find 
that temporal discretisation leads to softened vertex factors, modifies 
propagators and gives a natural cutoff for physical particle momenta. 



1 Introduction 

THIS paper is the third in a series devoted to the construction of discrete 
time classical and quantum mechanics, based on the notion that there is a 
fundamental interval of time, T. The objective is to investigate the proper- 
ties of a dynamics where continuity in time, and hence differentiability with 
respect to time, has been abolished. With no velocities, there are no La- 
grangians in the ordinary sense, and then there are no canonical conjugate 
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momenta or Hamiltonians either. It would appear then to be a catastrophic 
recipe for recasting the laws of classical and quantum physics, but as we 
try to show in this paper, this is not really the case. Moreover, there is 
every prospect for finding novel features of the dynamics not encountered in 
continuous time mechanics which may go some way towards alleviating the 
divergence problems encountered in conventional quantum field theory. 

The first paper of this series, referred to as Paper I fll], introduced basic 
principles for the temporal discretisation of continuous time classical and 
quantum particle mechanics. The second paper, referred to as Paper II 0, 
applied these principles to classical field theory, including gauge invariant 
electrodynamics and the Dirac field. These papers should be consulted for 
further explanation of our notation and methodology. In this paper, referred 
to as Paper III, we apply the techniques of Paper I to the quantisation 
of the scalar field systems studied in Paper II, i.e., we discuss quantised 
discrete time scalar field theory. 

Following the analysis of the earlier papers, we denote by V our process 
of discretising time using virtual paths and by Q the process of quantisation 
using transition amplitudes based on the system function, each of these pro- 
cesses being applied to some classical Lagrangian L. Then we can say that 
Paper I discusses models of type VL and QVL whereas Paper II discusses 
models of type DC, where £ is a Lagrange density. Now since such a density 
may be associated with the first quantisation of a classical theory, i.e., to 
QL models, as discussed in Paper II for the Schrodinger equation, models 
of type VC may be regarded as equivalent in some sense to those of type 
'DQL. This allows a direct comparison of the processes VQ and QD, and in 
Paper II it was argued that these are not the same in general. 

The present paper considers models of type QVC. Because such models 
may be regarded as equivalent in the above sense to those of type QVQL, 
then Paper III may be considered to be a discussion of discrete time second 
quantisation. Note however that the QVQ process used in this paper is not 
in general equivalent to the process VQQ because the V and Q processes 
do not commute. This means that our paper discusses the quantisation 
of discrete time classical field theories and not the temporal discretisation 
of quantum field theories, such as in lattice gauge theories. In the latter, 
discretisation is regarded as an approximation which becomes exact in the 
continuum limit. In our approach our mechanics is regarded as exact at each 
stage and the continuum limit is taken only to make comparisons with stan- 
dard formulations. This is a significant difference between our approach and 
various other formulations using a discrete time, because of our insistence on 
adhering to the principles of the formulation at all stages. In particular, the 
constants of the motion are constructed to be exact and not approximately 
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conserved up to some powers of T. 

An important question which arises naturally in the context of discrete 
time and/or space mechanics is that of Lorentz covariance. Our answer is 
that Lorentz symmetry of say scattering matrix elements emerges in the ap- 
propriate limit, such as T — >■ 0, and other than that, is not really something 
to worry about, as it is regarded here as an approximation to a deeper un- 
derlying structure. An analogy with representational art is useful here. If we 
liken continuous time theories to pictures drawn on normal canvas, then our 
discrete time mechanics is a picture drawn on a conventional analogue televi- 
sion screen. In the former model of spacetime it is frequently speculated that 
continuity might break down, perhaps at Planck scales (we do know that a 
real canvas is made up of atoms), but otherwise, continuity on the plane of 
the canvas exists at all levels and carries with it all the associated symmetries 
of the plane, such as translation and rotational invariance. On a television 
screen, however, we have two perspectives. From a distance, a television 
picture really does look hke one painted on a canvas, but a closer look would 
readily show the horizontal lines which make up the picture. There is a dis- 
creteness vertically, but a continuity horizontally. Likewise, in discrete time 
mechanics, there is a discreteness along the time axis with all the normal con- 
tinuity along the space axis. Like a television picture, there is less symmetry 
when viewed close up than when viewed at a distance, and it would be futile 
and in principle wrong to try to pretend that such long-distance symmetries 
should exist at all scales. What we are doing, therefore, is more like exploring 
the mechanics of a television set rather than the pictures drawn on it. This 
suggests that discretisation of time in the context of General Relativity is an 
obvious candidate for investigation. 

The art analogy can be pursued further. Discretisation of space as well 
as time, such as in lattice gauge theories and the work of authors such as 
Bender et al and Yamamoto et al @], gives a lattice space-time picture 
which corresponds to what occurs on a computer monitor, where the picture 
is fully digitised. This form of discrete space-time mechanics is inherently 
different to our discrete time mechanics and the two should not be confused. 

Another important question related to the issue of Lorentz covariance is: 
in which inertial frame are we discretising time? Of course, if we believed 
in an absolute time in the strict sense of Newton then we would have an 
immediate answer. However, we are approaching discrete time from a more 
modern perspective and the problem is a very real one for us. Our answer 
is to go beyond special relativity and consider cosmological perspectives. It 
is an empirical fact that, from the point of view of observers on the earth, 
we are moving at a speed of about 500 to 600 km per sec relative to a 
frame of reference in which the cosmic background radiation field of Penzias 
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and Wilson is isotropic P (the so-called dipole effect). According to the 
Cosmological Principle, we should be able to find a local inertial frame in 
the neighbourhood of each point in space time with the same property, i.e. 
one in which the cosmic background radiation field is isotropic to a very high 
degree, except for tiny ripples equivalent to those recently observed by COBE 
0. This frame should be unique at each point, up to spatial rotations. We 
will refer to this frame as the local absolute frame. 

If we allow that our universe is reasonably well described via a Robertson- 
Walker metric, then as we change position and time we expect the local 
absolute frame to change as well. However, it will always be empirically 
identifiable at each place and time in the universe. The standard co-ordinate 
time in such a frame is called the proper or co-moving time in the usual 
formulation of Robertson- Walker cosmology, and represents the proper time 
of a point particle (or galaxy) at rest relative to the local cosmic mass dis- 
tribution. In answer to the question posed above, we suggest that time is 
discretised via local absolute frames. 

Care should be taken to keep in mind that throughout this paper, when 
we discuss Minkowski spacetime and its temporal discretisation, we are really 
referring to local inertial frames. Of course, there is the additional question of 
local variations due to gravitational disturbances arising from locally inhomo- 
geneous matter densities. Answering this question amounts to constructing 
a discrete time analogue of general relativity, which will be reserved for a 
subsequent paper in this series. Since we are interested in applications to 
particle theory in this paper, we will not consider these issues further here, 
except to make a final observation about this line of thought. If we were dis- 
cretising space as well as time (which we are not), we would have to consider 
the additional question: how are we discretising space? If we were choosing 
the simplest sort of discretisation scheme, a cubic lattice (say), then we would 
have to specify three spatial orthogonal cartesian axes. Until recently there 
was no evidence of any spatial anisotropy on truly cosmological scales, so 
we had no criterion for picking out any special directions in space. We note 
however the very recent observation of the so-called corkscrew effect reported 
by Nodland and Ralston which if confirmed will require a re-assessment 
of the position. The comic background radiation field, however, does give 
us a working prescription for picking out a unique timelike direction at each 
point. 

Because of the relatively greater complexity of discrete time field theory 
compared with conventional field theory, we have restricted our attention in 
this paper to scalar fields. The general features found here should find their 
direct analogues with the Dirac and Maxwell fields. We reserve the further 
discussion of these fields to the next papers in this series. Our principal 
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aim in this paper is to discuss how the process of discretising time alters 
Feynman rules for scattering amplitudes and scattering cross-sections. Issues 
of renormalisation are left for later papers in the series. An important feature 
of the present investigation is the discrete time oscillator, which is directly 
related to free particle states used to define in and out states. 

In §2 we discuss the quantisation of scalar fields, using Schwinger's action 
principle to derive ground state expectation values of time ordered products. 
Then in §3 we apply these methods to the free neutral scalar field. The results 
are in agreement with the more direct calculation of the quantised discrete 
time harmonic oscillator discussed in Paper I. We examine in more detail 
the free scalar field propagator and the free field commutators, the results 
being consistent with the vacuum expectation values discussed previously. 
We find that the free particle creation and annihilation operators have a 
natural cutoff in physical particle state momenta, corresponding to what we 
call the elliptic regime. Although it is possible to construct linear invariants 
of the motion outside this regime, states associated with such operators have 
various expectation values which diverge or tend to zero in the infinite time 
limit, and this makes them unsuitable for representing physical particles. If 
T is the fundamental time interval and the discrete time analogue of energy 
E is defined by E = \/p.p + m'^ in natural units (where c = h = 1), then our 
formulation leads to the condition TE < -\/l2 for physical in or out particles, 
which we call the parabolic barrier. This barrier manifests itself in a number 
of ways. For example, the particle fiux density associated with each creation 
operator is found to be modified by a factor yJl-T^E^/12, which makes 
physical sense only in the elliptic regime. 

We then turn to interacting scalar fields theories. In §4 we set up the 
discrete time reduction formulae needed to calculate scattering cross sections 
and then discuss the pcrturbativc expansion of the vacuum expectation values 
of discrete time ordered products for a specific example, ip^ scalar field theory. 
We give the discrete time analogues of the Feynman rules in configuration 
space and in momentum space. In §5 we present a scattering calculation for 
the box diagram to illustrate the formalism and then give general rules for 
scattering amplitudes. Finally in §6 we give a number of applications of our 
scattering amplitude rules. We find that in each case there is a conserved 
quantity in scattering processes analogous to energy, related to the existence 
of what we call a Logan invariant of the system function. Fortunately, the 
LSZ formalism is powerful enough to reveal the existence of such a Logan 
invariant in a scattering process without the need for us to find it explicitly 
for the fully interacting system. 

Our analysis reveals that for ip^ interactions our discrete time Feynman 
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rules involve vertex softening in the basic diagrams, before any renormalisa- 
tion effects are considered. This may be a significant feature of more realistic 
interactions. Also, the propagators associated with internal lines are mod- 
ified and we use them to show how Lorentz covariance can emerge as an 
approximate symmetry of the mechanics. There is therefore some prospect 
of our programme making some progress towards the alleviation, if not com- 
plete removal, of divergences in the traditional renormalisation programme 
of continuous time relativistic quantum field theory. 



2 The discrete time quantised scalar field 

We turn now to the quantisation of the neutral scalar field. Following the 
methodology and notation discussed in Papers I and JJ, particularly the 
discussion in Paper I on the quantised inhomogeneous oscillator, the discrete 
time system function for a system with a scalar field (p degree of freedom 
coupled to a source j is chosen to be given by 

F"" [j] = F^+'-tJ rf^x + J , (1) 

where = / d^:Kj-'^ is the system function in the absence of the source. 
There are other ways of introducing sources into the system, but the above 
method was found to be most practical. Since these sources are eventually 
switched off, it does not really matter how they are introduced, as long as 
they are dealt with consistently according to the principles of discrete time 
mechanics. 

With the above system function the Cadzow equation of motion P, H is 

_l_{f» + f..-.} + r,„(x)^o, (2) 

which reduces to 

^ + - V . ^ + + n, ^ 0. (3) 

where JF" is the system function density in the absence of sources. Here we 
use the symbol = to denote an equality holding over a true or dynamical 
classical trajectory. 

The action sum A^'^^ [j] in the presence of sources for evolution between 
times MT and NT is 

N-1 „ 

+T j d'^Jn^n, M<N. (4) 
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Use of the discrete time Schwinger action principle |I| 

5{(l),N\%l),My =i{(t),N\5A^^ [iWilj.MY, M<N (5) 
leads to the functional derivatives 

■{<P,N\tij,My = U<p,N\0M{^)\^,My 



T SJm (x 

^jA^{<P,N\^,My = (0,iV|^„(x)|^,My, M<n<N 
T djn (x) 

(0,iV|^,MV = U<p,N\0^{^)\^,My. (6) 



T 5jN (x 
Also, we find 
—iS—iS 



T djn (x) T djm (y) 

+ {(j),N\0r. (X) 0m {y)\^,My 



(M <m,n < N) (7) 

where 9^ and 6n are the discrete time step function and discrete time delta 
defined in Paper I. We will write (|^ in the form 

^— — -(0, N\4j, My = {<p, N\f0^ (x) 0m (y) 1^, M)^ (M < m, n < iV), 
T (x) T djm (y) 

(8) 

where T denotes discrete time ordering as discussed in Paper I. 

In applications we will normally be interested in the scattering limit = 
—M — > oo and in matrix elements involving the in and out vacua. We shall 
restrict our calculations to such matters. This means we will discuss the 
r-point functions defined by 

, ^^ (0°"lt(^n, (Xi)...^„, (X,)|0-V- 

1 -i5 -iS 

Z [j] , (9) 



Z[j] T6jn, (xi)-T5j„, (x,) 
where 

Z [j] = (0°"*|0^'^)^' (10) 

is the ground state (vacuum) functional in the presence of the sources and T 
denotes discrete time ordering. 
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An important question here concerns the existence of the ground state. 
In common with most continuous time field theories, we have no general 
proof that a ground state exists for interacting discrete time field theories. 
Moreover, in discrete time mechanics there is no Hamiltonian as such, so 
the question becomes more acute. However, for free fields, there will be 
what we refer to as a compatible operator corresponding to some appropriate 
Logan invariant 0. This is the nearest analogue to the Hamiltonian in 
continuous time theory. Moreover, the appropriate compatible operator for 
free neutral scalar fields is positive definite and this allows a meaning for the 
in and out vacua to be given. 

3 The discrete time free scalar field 

3.1 The free scalar field propagator 

Given the continuous time Lagrange density 



Co = \d,^d^^ - i/iV' 




we use the virtual path approach discussed in to find the system function 
density 




In the presence of the sources we take 




(13) 



and then the Cadzow equation of motion is 




(14) 



We now define the spatial Fourier transforms 




(15) 
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and then the equation of motion becomes 

{ '^■-^:^""^ ^- '^""7^"'' }^.(p).j.(p). (16) 

where E = v^p^"j)+/? and f/„ is the classical temporal displacement operator 
defined by 

Unfn = fn+1 (17) 

for any variable indexed by n. The solution to (JT^ with Feynman scattering 
boundary conditions is 

oo 

(P) = (P) - 2^ E Ar'"(p)L(p), (18) 



m=— oo 



where Lp^^^ (p) is a solution to the homogeneous equation 

{ '^"-^r"'^ g' '^" + 7^"'> }^i°'(p)^0 ,19) 

and (p) is the discrete time Feynman propagator in momentum space 
satisfying the equation 

This equation for the propagator may be written in the form 

{Un - 2r]E + U-'}AI (p) = -TESn, (21) 

where 

= -TT—p^. (22) 



Using our experience with the discrete time harmonic oscillator propaga- 
tor discussed in Paper I, we may immediately write down the solution for 
the propagator in the form 

Af (p) = TTT^^e-^l"!^^ = {e-*"^-e„ + 5„ + c'^'^qA , (23) 

where rjE = cosBe- As discussed in Paper I, this expression holds for the 
elliptic and the hyperbolic regimes with suitable analytic continuation. In 
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the continuous time limit T — > 0, nT — > we recover the usual Feynman 
propagator in a spatially Fourier transformed form, viz; 

hm A^(p) = -±.{e-^^^e{t) + e^'^e{-t)} 

T->0, n^oo, nT=t Zll/ 



/ 



2E 



27r o;^ — p • p — //^ + ie 
= Kf (p, t)^ j d^:)c e-'P-'^Ap (x, t) . (24) 

Turning to the quantisation process, the functional derivative satisfies the 
rule 

^— ^Jm (y) = Sn-mS' (x - y) . (25) 

With the definition 

S 



we find 



and then 



Sjnip) (27r)'i 5jn(x) 

s 

Sjn (P) 



(26) 



Jm (q) = Sn-mS^ (p - q) (27) 



gives 



Hence we find 



(0-V„(x)|0-)^- = -^^-^ZL7l (28) 



(0-*|(^„ (p) |0-)^- = - '^y , f Z [j] . (29) 



Zo \j] = Zo[0] exp I -iiT^ f; / d'^d'yjn (x) Ar"^ (x - y) jm (y) | , 

I n,m=—oo J 

(30) 



where 

d^p 



(27r) 

This propagator satisfies the equation 

r - 2 + [/-I 



Af(x)= / e^P-A^ (p) . (31) 



2^2 



+ (^^ - V^) + ^-"^ I A- (x) = (x) ■ (32) 
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3.2 The free field commutators 



In this section we use (|30|) to obtain the vacuum expectation value of the 
free field commutators. Writing the propagator (|23|) in the form 



(p) = CE [^-"e„ + 5n + z^Q 
then we find for the elliptic and hyperbolic regimes 



(33) 



ce 



EV12-T^E^' 
-3 



TE < V12 
TE > VU. 



EVWW^U' 
An apphcation of (|^) gives 

{O\0n+i (x) 0n (y) |0) = tA], (x - y) , 
from which we deduce 

^3p eip{x-y) 



(34) 



(35) 



(0|[^„,+i (x),<^„ (y)]|0) 



'6iT 
-6i 



4:7iT |x — y| 



(27r)^ 6 + (p ■ p + /i2) T2 

g-^M2+6/T2|x-y|_ 



(36) 



Both elliptic and hyperbolic regions of momentum space contribute to this 
result, which has the form of a Yukawa potential function. 

We turn now to the direct approach to quantisation discussed in Paper I. 
If we define the momentum 7r„ (x) conjugate to (pn (x) via the rule vr^ (x) = 
^ rF" then (pi) gives 



<5l/3„(x) " 



TTn = + ^ ^ ) (^"+1 + 2*^^ 



(37) 



for the free field. The naive canonical quantisation discussed in Paper I is 
equivalent in field theory terms to 



from which we deduce 

[0n+l (x) , 0n (y)] 



[7r„ (x) , 0n (y)] = -iS^ (x - y) , 

d^p gip-(x-y) 



(38) 



-6iT 
—6i 



47rT |x — y 
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(27r)^ 6 + (/i2 + p ■ p) T2 

g-V/^2+6/T^|x-y|^ (39) 



assuming [ipn (x) , (p^ (y)] = 0. This is consistent with the approach to quan- 
tisation via the Schwinger action principle from which we obtained (|36D . 

We note that the reason this works is that the system function for a free 
field is an example of what we call a normal system in Paper I. For interact- 
ing field theories this will no longer be the case and then the commutators 
analogous to the above will probably no longer be c-functions. We recall 
that in continuous time field theories, interacting field commutators are not 
canonical in general either, so the analogies between discrete and continuous 
time mechanics hold well here also. 

For the free field particle creation and annihilation operators we define 
the variables 



an (p) = zF^^e^"^- [(pn+i (p) - e''^(pn (p)J , (40) 

where 

6T 



fiT r 
Q^T2E2 ^ (P) = j rf^xe-P>„ (x) (41) 



and the momentum p is restricted to the elliptic region. Then we find 

a„ (p) , at (q)] = 2E^l - T^E^ /12 {2^f (p - q) (42) 

when we quantise and use (PU]). This tends to the correct continuous time 
limit as T — i> 0. 



If we interpret the factor 2Eyl — T^E'^/12 in the above as a particle 
flux density then this will be indistinguishable from the conventional density 
2E in continuous time field theory for normal momenta, but falls to zero 
as the parabolic barrier TE = \/l2 is approached from below. This sug- 
gests that there is in principle a physical limit to the possibility of creating 
extremely high momentum particle states in the laboratory or of observing 
such particles in cosmic rays. This should have an effect on all discussions 
involving momentum space, such as particle decay lifetime and cross-section 
calculations, and in the long term, on unified field theories. 



4 Interacting Discrete Time Scalar Fields 
4.1 Reduction formulae 

In applications to particle scattering theory we shall be interested in incom- 
ing and outgoing physical particle states, with individual particle energies 
satisfying the elliptic condition TE < \/T2. We note that energy is defined 
here via the linear momentum p by the rule 

E = +V^MHV- (43) 
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This is the only meaning we give to the term energy. 
Given the annihilation and creation operators 



(fn+l (x) - e*^^(^„ (x 



a„(p) = ir^^y (i3xe'"^^-'P-^ 
a+(p) = -irii|d^xe-^"^-+^P-^[(^„+i(x)-e-^^-(^„(x)] (44) 



where 



6T 



cos 9E = r]E^ 



6 - 2r2£;2 



(45) 



6 + r2£;2 ' ^ - 6 + T2£;2 ' 

then a direct apphcation of the standard reduction formahsm gives the re- 
duced matrix elements 



oo „ 

^ J d3xe-^'^^-+^P-^i?;::;(a°"*|TC^„ (x) 



(46) 



and 

(a-*|a,„,(p)fC|/?-)R 



z ^ /d=^xe^"^--^P-^:^(a°«*|T'C(^„(x)|/?-), 

n=— oo 

(47) 

where C denotes any collection of field operators and 

^ = - 2riE + t/-') . (48) 

Using these results we can readily write down the scattering matrix for a 
process consisting of r incoming physical momentum particles with momenta 
Pi, P2, Pr and s out- going particles with momenta qi, q2, q^. 



4.2 Interacting fields: scalar field theory 

We turn now to interacting scalar field theories based on continuous time 
Lagrange densities of the form 

C^Co-Viif). (49) 

In order to illustrate what happens in discrete time quantum field theory, we 
shall discuss the details of a scalar field with a ip^ interaction term, deriving 
the analogue of the Feynman rules. 
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In the presence of sources the above Lagrange density leads to the system 
function 

^" [j] = ^(0) - T d\ J t/^X V i^n) + \T j d^^ + Jn+lVn+l} , (50) 

where we use the virtual paths 

ifn (x) = U^ipr, (x) = Xipn+i (x) + Av5„ (x) , < A < 1, A = 1 - A, (51) 

as discussed in Paper II for neutral scalar fields. Here and below we shall 
find the operator 

f/„^ = XUn + X (52) 

particularly useful, where [/„ is the classical temporal displacement operator 
defined by (|l^. The vacuum functional is now defined via the discrete time 
path integral 

Z[j] = I [d^]exp{tA[j]} 

= f[ if [d^n])exp{zA[j]}, (53) 

n=-oo ^-^^ ^ 

where 

oo oo „ 

Mj]^ E ^'^b]= E Fio)[j]-4^ Vi^n) (54) 

n=-oo n=-oo JnXx 

and the (/?nare functionally integrated over their spatially-indexed degrees of 
freedom. In the above and subsequently we shall use the notation 

i =T f^dX [ d^x (55) 

JnAx Jo n=-oo-^ 

whenever such a particular combination of spatial integration, summation, 
and virtual path integration occurs. This replaces the four- dimensional inte- 
gral / d'^x = J dtd^x found in normal relativistic field theory. 

We now postulate our quantum dynamics to be governed by the equation 



exp{!/l|j]} = 0, (56) 



which is equivalent to a vacuum expectation value of the Heisenberg operator 
equations of motion derived formally from Cadzow's equation (^ . Integrat- 
ing by parts, we arrive at the more convenient expression 

Z [j] = exp I -if^^^ V (Vnx^) I Zo [j] , (57) 
14 



where 



{oo „ 
n,m=—oo 



mKy) 

(58) 



and 



T "5i„(x) 

= -|a^-^ + A^|. (59) 

Turning now to (fHheory, we recall that with hindsight the potential 
((^) is normally taken to have the form 

v^'^ {<P) = I - r^} , (60) 

where the (infinite) subtraction constant F is formally given by 

r = 3iAF(0). (61) 

This has the role of cancelling off self-interaction loops at vertices in the 
Feynman rules expansion programme. We find that for discrete time, the 
same effect is achieved by taking the potential to have the form 

V^'^ i^n) = I {<^n - r^^n} (62) 

where 

f = 2iA% (0) + y [a]. (0) + A;^^ (0)] . (63) 

The first objective is to find a perturbative expansion for Z [j] , which we 
write in the form 

Z\j]^Zo\j]+Z,\j]+Z,\j] + ... (64) 

where 

Z,\j] = --i V^'^VnX.)Z,^i\j], p=l,2,... (65) 

Having found Z [j] we then calculate the required vacuum expectation value 
of time ordered products of fields by functional differentiation in the standard 
way. The results lead to a set of rules for a diagrammatic expansion analogous 
to the Feynman rules in continuous time theory, with specific differences. The 
details of the calculations are omitted here as they are routine and tedious, 
but the results are as follows. 
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4.3 Feynman rules for discrete time-ordered products 

The objective in this subsection is to present the rules for a diagrammatic 
expansion of scattering amphtudes in the absence of external sources. The 
latter are used merely to provide an internal handle on the correlation func- 
tions of the theory and are set to zero at the end of the day. This programme 
is carried out in two stages. In this subsection we give the rules for the eval- 
uation of successive terms in a Feynman diagram type of expansion for the 
vacuum expectation value of the time ordered product 

(0°"*|f ^1 (xi) <^2 (X2) ...0k (xfc) |0*") (66) 

with k discrete time scalar fields; we shall give the rules for a system with 
interaction given by (p^), so the expansion is effectively in powers in the 
coupling constant g: 

1. first find the ordinary continuous time Feynman rules in space-time; 

2. draw all the different diagrams normally discussed in this programme; 

3. for a given diagram with V vertices and I internal lines find its con- 
ventional weighting factor u, such as the well-known factor of ^ for the 
simple loop in ip^ theory; 

4. at each vertex, associate a factor 

tgri =igT['d\ f; /d^z; (67) 

JmXx Jo m=-oo 



5. for each external line running from the external point (n, x) to a vertex 
with indices (m, A, z) assign a propagator 

^t/^A™-" (z - x) ; (68) 

6. for each internal line running from vertex (mi, Ai, Zi) to vertex (m2, A2, Z2) 
assign a propagator 

^^^\f/™^.AF™(z2-zO; (69) 

7. do the A integrals. 
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It is in general much more convenient to perform the virtual path integra- 
tions (over the A's) after the diagrams have been written down rather than 
before the diagrammatic expansion. In many cases the operator acting 
on an external propagator can be transferred to act on internal propagators 
using the rule 



E [U::^fm)9m= E /-f^^^?-, (70) 
m=— oo m=— oo 

for any indexed functions /„, (yf^, where we define 

U^g-m = >^U;^^grn + ^9m = Hm-l + ^9m- (71) 

However, this does not work so conveniently whenever two or more external 
lines meet at the same vertex. 

To illustrate these rules in operation consider the conventional pertur- 
bation theory expansion of the time ordered product (0|T(^ (xi) (p (^2) |0) in 
powers of the coupling constant. The conventional Feynman rules give the 
expansion 

(0|T(^(X1)(^(X2)|0) = i/\f{Xi-X2) 

-\g^ j d'^zid'^Z2Ap {xi - Zl) Ap {zi - Z2) x 
Ap {Z2 - Zl) Ap {Z2 -X2) + (g^) . (72) 



The second term on the right hand side corresponds to the single loop dia- 
gram with V = 2, I = 2 m. ip^ scalar theory and is divergent. Part of the 
motivation for investigating discrete time field theory is the hope that the 
corresponding diagram might be modified in some significant way. 

Using the rules outlined above the analogue expansion in discrete time 
gives 

(0|T<^„,(xi)<^„,(x2)|0) = 2A^^-"Mxi-X2) 

-\H f {f/^\A---(zi-xi)}x 

v/mi AlZi v/m2 A2Z2 

{ui\U^^A-^-^ (Z2-Z1)} {u^,Ui\AT--^ (Z1-Z2)} 
[ut\A^^--- {^2-^2)} + O {g') . (73) 

For this particular process the second term on the right hand side can be 
rewritten using the rule ([70|) to give 

(0|f(^„,(xi)<^„,(x2)|0) = zA^i-'^Mxi-X2) 



X 



-\g'± i A^-"^ (zi-xi) X 

Jmi Aizi J m2A2Z2 
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^ mi ^ 1712 



TtM TtM /\m2-mi / ^ 



X 



A^-"^ (Z2-X2) + O (g') , (74) 



using the symmetry 



A^ (x) = A^" (-X) . (75) 



The integrals over fi and A can be integrated at this stage to give a mul- 
titude of subdiagrams distinguished by different split times, which is the 
ultimate effect of the discretisation process. The various subdiagrams con- 
tributing to the loop diagram are shown in Figure 1, each with a numerical 
factor. The sum over all numerical factors for this diagram should add up 
to 144. The full amplitude corresponding to the loop diagram is the sum of 
each of these sub-diagrams, times the numerical factor for each sub-diagram, 
divided by 288, taking into account the original weighting factor of one half. 
By using symmetry arguments it can be shown that the twenty nine distinct 
diagrams in Figure 1 reduce to the twelve diagrams shown in Figure 2. 

The above rules are relevant to vacuum expectation values of discrete 
time-ordered products of field operators. For particle scattering matrix ele- 
ments the rules become simpler, as discussed next. 



5 Scattering amplitudes 

We are now in a position to discuss particle scattering amplitudes. First 
we explain how the scattering amplitude for a two-two particle scattering 
process based on the box diagram. Figure 3, is calculated, and then we shall 
state the results for the general scattering diagram. This diagram was chosen 
because it involves a loop integration. 

5.1 The two- two box scattering diagram 

Consider two incoming scalar particles with 3— momenta a, b respectively 
scattering via a the box scattering diagram shown in Figure 3, into two outgo- 
ing particles with 3— momenta c and d respectively. Each of these particles is 
associated with a 9 parameter as given by (|^) which lies in the physical par- 
ticle interval [0, vr). Negative values of such a parameter correspond to waves 
moving backwards in discrete time and would be interpreted in the usual way 
as anti-particles in the Feynman-Stueckelberg interpretation. Both positive 
and negative values occur in the discrete time Feynman propagators, just as 
in conventional field theory. 
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Using the reduction formulae in §4.1 we may write for the scattering 
reaction amphtude Sif 

% ^ (0''"*|a„„i(d)a„„,(c)a+ (b)a+ (a)|0'")K 



„-4 



4 oo . 

n„5: / 



d^x,- e 



— mi^a+ia-xi — m2Sb+ib-X2„m36c— ic-X3 in46d—i<i-X4 



X 



i j=l nj=— oo ' 
->-rT >rrT >rrT > 



-f^'ni,a-f^n2,b-f^n3,c-f^n4,d(0°" \T(fini (xi) (^^^ (xs) (^„3 (xg) (^^^ (X4) |0'"), 



where 

with 
r(a) 



6T 



kZ!. = (a) - 277 (a) + C/-^) 

6 - 2T^E'^ 



6 + T^E^ 



a-a + /x^, ?7(aj = 



6 + T'^El 



(76) 
(77) 

COS^a, (78) 



and similarly for the other particles. 

Next we expand the 4-point function according to the rules outlined in 
§4.3 and consider for the purposes of this discussion only the contribution 
associated with the box diagram of Figure 3, viz 



(O^^'lTf^n, (Xi) ipn, (X2) Vn,(^^)^nM^W)\BOX 

Ui^Ui^I^T"^^ (Z4 - Z3)l teC/^lA--- (Zi - Z4) 



mz-m2 I \ 



X 



(79) 

The next step is to do the Xj integrals, converting the two point functions 
on each external leg of the diagram to its momentum space form, using 



An 



(p) = y dxe^P-^A^ (x). 



(80) 



Then we use the result 



K„,pA^ (p) = -K 



taking care to bring the operators and summations into the brackets when- 
ever the operators occur. This effectively amputates the external legs of 
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the diagram. Then we can immediately carry out the summations over the 
external integers rij and arrive at the simplified form 



S^f = {9T)'(y[ E [d\J 



(i^z,- e 



ia-zi+ib-Z2— ic-Z3— id-Z4 



X 



(82) 



{U^^^e-^^^o.} {U^\e-^^''-] {f/„^;e-'3e'=} {u^\e^^^'-} x 
UtUt\/^7'--' (Z2-Z1)] [C/A^3C/^^,A--- (Z3 - z^)] X 

ui^u^JAT'"^' (^4 - Z3)] [u^ui^^^-^^ (zi - Z4); 

Now we use the representation of the propagator 

(x) = ^ / ci^k £ ^e-'"^+^''-A^(k, ^) (83) 
and evaluate the z^integrals to find 

% = /(27r)^5-^(a + b-c-d) (n E fd\~rdeAj 

I^Mg-^mi^a} {c/A.e--^^'^} {C/A|e^"^^^=} {c/^:e^-^^<^} x 
[f/^^.^™\e-("^)^^ A^ (k,^i)] [f/^33f/^^^e-(-3-.) (k + b, 9,) 



(27r)^ 



f/™lf/;^i:.e-^(™^-™3)A^ (k + b - c,^3) 



X 



(84) 



Here we see the appearance of overall linear momentum conservation, as 
expected. Next we use the result 



tt\ imd „imd 



where 
to find 



(85) 
(86) 



/ 4 00 „i 1 r-^ \ 
S,f = g'{2nfS'{^ + h-c-d)i]l E J, ^^^^ J_J^A 

yj—lrrij—co J 



X 



t (7^k 



20 



X 



^i(m,-ms)e,f^^ (^2) /a*3 (^2) (k + b, 62) X 

^i(ms-m,)e,f^^ (^3) (^^3) A^ (k + b - C,^3) X 
^iim,-m,)e,f^^ (^4) /m (^4) A^ (k - a,^4) . 

(87) 

We arc now able to do the summations over the mj. We notice that each 
summation gives a Fourier series representation of the periodic Dirac delta, 
viz 

00 00 

J2 e^"'^ = 27r J2 S{x + 2rmT) = 2nSp{x), 

m=— 00 m=—oo 

and so we find 

i2nr 

fx, (9) fl (9) ik,9) fx, {9 + 9,) f^ {9 + 9,) A^ (k + b, ^ + 9,) x 
/a3 {9 + 9,-9,)fl{9 + 9,-9,)Kp{\^ + h-c,9 + 9,-9,) x 
fx, {0 - 9a) fl {9 - 9a) A^ (k - - 9^) . 

(89) 

The crucial significance of this step is that we see the appearance of a con- 
servation rule for the parameters 9. This is despite the non-existence of a 
Hamiltonian in our formulation and the fact that we have not constructed a 
Logan invariant for the fully interacting system. 

We may go further and do the Ajintegrals. We define the vertex function 

V{9a,9,) = /'dA/*(^a)/A*(^b)/A(^a + ^b) 

Jo 

COS (6'a + 6'b) + COS (6'a) + COS (6'b) + 3 

- ^ (90) 

and so get the final result 

% = g'{2n)Up{ea + 9,-9,-9d)5^{a + h-c-d)x 

/(j^ k /■"" 
^ J J9V {9a, -9) V {9,, 9)V{9 + 9,, -9,) V {-9,, 9a- 9) x 

Kp (k,^) (k + b, ^ + 9^) A^ (k + b - c,9 + ^6 - 9,) Kp (k - a,^ - 9a) . 

(91) 
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A diagrammatic representation of the above shows that ^-conservation occurs 
at every vertex. 



5.2 The vertex functions 

The vertex functions V(0i^02) represent a degree of softening at each vertex 
arising from our temporal point sphtting via the system function. At each 
vertex the sum of the incoming 9 parameters is always zero, including inside 
loops, so the vertex function always depends on two parameters only. If 
we had a y?^ interaction we expect the vertex function will depend on three 
parameters, and so on. A graphical presentation of the vertex function is 
given in Figure 4. The vertex function has a minimum value of one quarter 
and attains its maximum value of unity when the 9 parameters are each zero. 
This corresponds to the continuous time limit T ^ 0. 



5.3 The propagators 

The propagators used in the final amplitude (|9TD are readily found using the 
basic definition 

oo 

(p, 9)^T e*"'^^ (P) (92) 

?i=— oo 

and the equation 

V„ - 2r]E + U-^) (p) = -r^;^,. (93) 



Then 

2(cos^-r/i5)A^(p,^) = -Tr^. (94) 

Now we need to choose the correct solution for Feynman scattering boundary 
conditions. This is done by referring to the Feynman —ie prescription, which 
corresponds to the replacement of E"^ in the above by E"^ — ie. This in turn 
corresponds to the replacement 

Ve ^Ve + ie. (95) 
Hence we arrive at the desired solution 

AAp,0) = ^. n ^ -y (96) 
/ (cose' — rjE — te) 

which holds for both the elliptic region — 1 < r^^; < 1 and for the hyperbolic 
region —2 < tje < —1. It may be verified that the indexed propagators (p3D 
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are given by the integrals 



-ind 



Af(p) = dee- 

j-e ^ 

rio' 

(97) 



2 [cosO — rjE — ie) 
Te P dz 



2m J _ 2(^r]E + ie)z + 1) 

the contour of integration being the unit circle in the anticlockwise sense. 
We find for example 

A" (p) = ^""^ e-^l"''^ (98) 
^^^^ 2ismeE 

in the case of the elliptic regime, T'^E'^ < 12, and 

(—1 "1"+^ r 

(P) = \ \ "" e-H^^ (99) 
2 smh7£; 

in the hyperbolic regime, T'^E'^ > 12. Here we make the parametrisation 

6 — 2T^ E'^ 

cos iC) ^r^E = y,^, , (100) 

where is a complex parameter running just below the real axis from the 
origin to vr (when ( is written as 9e) and then from vr to tt — iln(2 + -\/3) 
(when ( is written in the form vr — i'jE)- 

If in (|96|) we introduce the variable po related to 6 by the rule 

6 - 2p'^T'^ 

^^^^ = ~^, — ¥Pr^ sign (9) = sign (po) , (101) 

then we find 

Af (P,^) = / + 2 2^- V ^102) 

Po — P ''T^ + [Pq — — + «ej 

an exact result. From this we see the emergence of Lorentz symmetry as 
an approximate symmetry of the mechanics. If po in the above is taken 
to represent the zeroth component of a four-vector, with the components 
of p representing the remaining components, then we readily see that the 
first term on the right-hand side of ( |102| ) is Lorentz invariant. The second 
term is not Lorentz invariant, but we note it is proportional to T^. If, as we 
expect, T represents an extremely small scale, such as the Planck time or 
less, then it is clear that Lorentz symmetry should emerge as an extremely 
good approximate symmetry of our mechanics. 
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5.4 Comments 



The significance of our results is tliat not only is spatial momentum conserved 
during a scattering process, as expected from the Maeda-Noether 0, || theo- 
rem , but the sum of the 9 parameters of the incoming particles is conserved. 
This is the discrete time analogue of energy conservation, since in the limit 
T — > we note 

Or 



^m^ = E^ = ^Y>■V + ^^'■ (103) 

The 9 conservation rule is unexpected at first sight in that we have not 
discussed as yet any Logan invariant for the full interacting system function. 
It appears that the analogue of energy conservation occurs here because of the 
way in which we have set up our incoming and outgoing states and allowed 
the scattering process to take place over infinite time. The result would 
probably not hold for scattering over a finite time intervals, which would be 
the analogue of the time-energy uncertainty relation in conventional quantum 
theory. In essence, the LSZ scattering postulates relate the Logan invariant 
for in-states to the Logan invariant for the ont-states in such a way that 
knowledge of the Logan invariant for the intermediate time appears not to 
be required. This is true of the scattering formalism as we have demonstrated, 
but the bound state question would be a different matter. 

Although the conservation of 6'— parameters during scattering processes 
comes as a surprise it is a welcome one. Before the calculations were done 
explicitly, it was believed that the energy conservation rule in continuous time 
scattering processes would only arise in the limit T ^ 0. Such a phenomenon 
was discussed by Lee fllj in his discrete time mechanics, which differs from 



ours in that his time intervals are determined by the dynamics. That there 
is an exact conservation rule for something in our discrete time scattering 
processes regardless of the magnitude of T is an indicator of the existence of 
some Logan invariant. The surprise is that the something turns out to be the 
sum of the incoming 9 parameters, which suggests that our parametrisation 
of the harmonic oscillator discussed in Paper I was a fortuitously good one. 

We point out here that our parameter 9 is really an angle, unlike con- 
ventional energy or po in the above, and there is an implied periodicity. 
However, because there is no concept of Hamiltonian or energy in our theory, 
this periodicity does not cause any physically relevant side effects. Incoming 
or outgoing particles will be on-shell in the sense that their associated 9 pa- 
rameter can be restricted to takes a value in the interval [0,7r). Given this. 



then we may invert ( |101|) to find 

9 = p,T-^pIt' + 0{t''). (104) 
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Another welcome feature is the modification of the propagators and the 
appearance of vertex softening in the scattering diagrams. A detailed dis- 
cussion of the effects these features have on the divergences of various loop 
integrals found in the conventional Feynman diagram programme will be 
reserved for a subsequent paper. 

The scattering amplitude found above for Figure 3 reduces to the correct 
continuous time amplitude in the limit T — > 0. 

5.5 Rules for scattering amplitudes 

We are now in a position to use our experience with the box diagram Figure 3 
to write down the general rules for scattering diagrams. Consider a scattering 
process with a incoming particles with momenta pi,p2,...,Pa respectively, 
and b outgoing particles with momenta qi,q2,---q6 respectively. Make a 
diagrammatic expansion in the traditional manner of Feynman. For each 
diagram do the following: 

1. at each vertex, conserve linear momentum and 6 parameters, i.e., the 
algebraic sum of incoming momenta is zero and the algebraic sum of 
the incoming 9 parameters is zero; 

2. at each vertex associate a factor 

igTV{e,,e2), (105) 
where 9i and 6*2 are any two of the three incoming 9 parameters; 

3. for each internal line carrying momentum k and 9 parameter, associate 
a factor 

iT-'Ap {k,9) ; (106) 

4. for each loop integral, a factor 




5. an overall momentum-^ parameter conservation factor 

(27r)^ 6p {9p, + ... + 9pa- 9q, - ... - 9gb) 6^ (pi + ... + p„ - qi - ... - q^) 

(108) 

6. a weight factor lv for each diagram, exactly as for the standard Feynman 
rules. 
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6 Examples 



We are now in a position to give a number of examples of scattering amplitude 
calculations using the above rules. We restrict our attention to ip^ theory as 
an illustrative example. QED and the associated discrete time Feynman rules 
will be the subject of Paper VI in this series. 

6.1 Figure 5a: 

Consider the basic single vertex diagram of figure 5a with particles with 
linear momentum a. b fusing to form a particle with linear momentum c. 
Overlooking the fact that this process gives zero for on-shell momenta our 
discrete time Feynman rules give 

S^a = igT (27r)Up(ea + e,-e,) 5^ (ai + h-c)V(ai,h). 

(109) 

6.2 Figure 5b 

This diagram has a single loop. We find 

V {Oa, -&) V {Ob, e-ec)v {-e^, e) Af (k, e) (k - c, ^ - e^) Af (k - ^,9 - Oa) . 

(110) 

The nature of this diagram will be discussed in detail in subsequent papers. 

6.3 Figures 5c,d,e 

The order g'^ two-two scattering diagrams figures 5c,d,e give 

-^scde = -g^T{2Tif 5p{ea + Ok- e,-ed)5\si + h- c-d)x 
{V {Oa, -Oc) V {Oh, -Od) Ap (a - c,ea - 9,) + 
V [OaA) V i-e,, -Oa) (a + hA + 0^) + 
V {9a, -9a) V {9b, -9,) Ap (a - dA - 9^)} . 

(Ill) 
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6.4 Figure 5f 

This diagram is an example of a higher order tree diagram process involving 
no loops. We find 

S^f = ig^T{27r)Up{ea + eb-e,-ea-e,)6'{ai + h-c-d-e)x 

V {Oa, -e,) V {Ob, -6d) V (Oa-e^A-Od) x 

(a -c,ea- e,) (d -h,ea- e,) . (112) 

6.5 Figure 5g 

This diagram has a simple propagator loop and gives 

Af (a + h,ea + 9,) I / (i^k /' dO V {Oa, 9,) V {9^ + 9,, -9) x 

V {9, -9a - 9i,) V {-9c - 9d) Af (k, ^) A (k - a - b, ^ - - 9^)} x 
A^(a + b,^„ + ^fe). (113) 

The question of the divergence of this integral will be reserved for a subse- 
quent paper. 

7 Concluding remarks 

The application of the principles outlines in Papers I and // to scalar field 
theory has indicated that the conventional programme of constructing Feyn- 
man rules for scattering amplitudes goes over well into discrete time. Of 
course there are differences, and it is to be hoped that some of these will al- 
leviate if not overcome some of the divergence problems of the conventional 
field theory programme. An important point is that there occurs in our ap- 
proach a natural scale provided by T. It is possible that this will provide a 
renormalisation cutoff scale which will not have to be introduced by hand. 
Issues of renormalisation and divergence will be discussed in a later paper. 

A particularly important result which was not anticipated before the di- 
agrams were calculated is the conservation of the total 9 parameters over a 
scattering process. This occurs even though no Logan invariant correspond- 
ing to the total Hamiltonian has been found for the fully interacting theory. 

An important point to consider is the question of relativistic covariance. 
Clearly our process of temporal discretisation breaks Lorentz covariance, and 
with it the Poincare algebra. However, it should be admitted by any critic 
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that there is actually no empirical evidence that Special Relativity holds all 
the way up to infinite momentum. It is only an abstraction from limited 
experience that it does. Therefore, the Poincarc algebra has no more than 
the status of a really useful synthesis of limited experience. By requiring our 
parameter T to be small enough we should be able to reproduce all of the good 
predictions of continuous time mechanics, with the possibility of alleviating, 
if not removing, those aspects which are known to cause problems, such as 
divergences in the renormalisation programme. Moreover, we have given a 
principle based on the cosmic background radiation field for finding a unique 
local inertial frame in which time is discretised. 

Finally, if our discrete time programme could be caught out in a fatal way, 
then we would have what amounts to a proof that time is really continuous. 
This in itself makes our investigation a worthwhile one. 
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Figure 1: An expansion of the single loop diagram in configuration space. There 
is a total of 144 contributions, many of which can be equated and so we give the 
weighting for each type of contribution. The apparently disconnected pieces in most 
of these types arise because we have effectively split points in time when we wrote 
down the system function. In the limit T — these splits would collapse and all 
of these diagrams then reduce to 144 copies of the first type, corresponding to the 
conventional single loop diagram. 
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Figure 2: By using symmetry arguments, the various kinds of contributions shown 
in Figure 1 can be simphfied to the 12 varieties shown here. Again, the sum over all 
weightings is 144. 
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Figure 3: The two-two box scattering diagram used to illustrate the construction 
of the discrete time diagrammatic rules. 
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Figure 4: A plot of the vertex function V{x,y), where each argument hes in 
the interval [— 7r,7r]. The vertex function takes its maximum value of unity when 
X = y = 0, which corresponds to the continuous time limit T — 0. The vertex 
function never falls below the value 1/4. 
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Figure 5: A number of diagrams to which our discrete time scattering rules have 
been apphed. 
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